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We investigate high-order harmonic spectra from aligned diatomic molecules in intense driving 
fields whose components have orthogonal polarizations. We focus on how the driving-field ellipticity 
influences structural interference patterns in a macroscopic medium. In a previous publication [Phys. 

Rev. A 88, 023404 (2013)] we have shown that the non-vanishing ellipticity introduces an effective 
dynamic shift in the angle for which the two-center interference maxima and minima occur, with 
regard to the existing condition for linearly polarized fields. In this work we show through simulation 
that it is still possible to observe this shift in harmonic spectra that have undergone macroscopic 
propagation, and discuss the parameter range for doing so. These features are investigated for H 2 
in a bichromatic field composed of two orthogonally polarized waves. The shift is visible both in 
the near- and in the far-field regime, so that, in principle, it can be observed in experiments. 
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I. INTRODUCTION 


Recently, driving fields with nonvanishing ellipticity 
have attracted a great deal of attention as potential at- 
tosecond imaging tools M- In particular, with re¬ 
gard to molecules, these types of fields provide access 
to a whole parameter range which would not be ac¬ 
cessible otherwise. For instance, they may be used to 
probe degenerate orbitals, molecules that are difficult to 
align, and also allow the reconstruction of molecular or¬ 
bitals from a single-shot measurement. They can also 
be employed in molecular imaging of randomly oriented 
molecules II- 

This is made possible due to the physical mechanism 
behind high-order harmonic generation (HHG), namely 
the laser-induced recombination of an electron with its 
parent ion 10- By using driving fields composed of or¬ 
thogonally polarized waves, one can in principle control 
the angle with which the electron leaves, or subsequently 
returns to the ion [J, 2;, Q . This control has been real¬ 
ized experimentally by employing fields with frequencies 
io and 2 to and changing the relative phase between both 
waves H. 

A particularly interesting example is related to HHG 
in aligned diatomic molecules. Even if HHG is mod¬ 
eled within a very simple framework, namely the strong- 
held approximation (SFA) and the single-active electron 
(SAE), single-active orbital approximation, driving fields 
with non-vanishing ellipticity lead to interesting features. 
One of such features is the blurring and the splitting that 
occur in the two-center interference minima, which has 

been reported in Bairn. For linearly polarized fields, 

these minima are purely structural and have been related 
to electron recollision at different centers in the molecule 


(for the first discussion of this minimum see 0, [Til an d 
for reviews see, e.g., [lTj and our recent publication fl5|). 
If, in contrast, the field is composed of orthogonally po¬ 
larized waves, the interference condition changes. Indeed, 
in previous work, we have derived a modified two-center 
interference condition which incorporates the electron’s 
angle of return as a dynamic shift 11]. This shift is crit¬ 
ically influenced by the orbit along which the electron 
returns, and by the harmonic energy. Hence, a coherent 
superposition will lead to the above-mentioned blurring 
and, in extreme cases, to the splitting. 


A legitimate question is, however, whether this shift 
can be observed under more realistic conditions, such as 
experiments. In 11] we have anticipated a series of diffi¬ 
culties that may arise. First, even at the single-molecule 
response level, there are coherent superpositions of or¬ 
bits, which will be detrimental for the identification of 
individual shifts. Second, the core dynamics may mask 
these shifts further and has not been incorporated in this 
model 0 , 0 . Third, these shifts depend on the driving- 
field intensity, which will vary across the beam profile. 
Fourth, there is no evidence that these shifts survive ef¬ 
fects that occur during HHG propagation in a macro¬ 
scopic medium. 


On the other hand, propagation may facilitate the ob¬ 
servation of these shifts, since the long and the short or¬ 
bits in the dominant pair phase match differently 00 - 
Thus, by carefully choosing the propagation conditions 
one may be able to select the contributions, and the 
shifts, associated with individual orbits. This type of 
selection has been successfully used in the past to op¬ 
timize attosecond-pulse production (for a seminal pa¬ 
per and a comprehensive review on this topic see, e.g. 
[18i |. and 0)- In fact, it has even been shown that to 
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obtain clear enough interference patterns between both 
orbits requires a very delicate tuning on the propaga¬ 
tion conditions [U [22]. Further trajectory selection 
may be achieved by creating temporal or spatial gates, 
by playing around with the pulse shape, polarization 
and/or macroscopic conditions. Examples of different 
trajectory-selection mechanisms are provided in the re¬ 
views gain. 

In this paper, we perform a detailed analysis of 
how these shifts manifest themselves in a macroscopic 
medium, for individual orbits. We verify that the shifts 
behave in different ways for the long and short orbits, and 
are relatively stable across the beam. This implies that 
they are not washed out by propagation effects, but, in 
principle, can be seen in a more realistic, macroscopic sce¬ 
nario. In our investigations we consider H 2 to avoid mul¬ 
tielectron effects associated with the core dynamics. We 
consider both linearly polarized fields and bichromatic or¬ 
thogonally polarized fields of frequencies ui and 2 to. For 
the single-molecule response, we employ the strong-field 
approximation (SFA) and the steepest descent method. 

This paper is organized as follows. In Sec. El we bring 
the necessary background to facilitate the subsequent dis¬ 
cussion. The results are presented in Sec. uni from the 
overall behavior of the effective shifts with the driving- 
field intensity to an analysis of how they behave nracro- 
scopically. Finally, the main conclusions to be drawn 
from this paper are stated in Sec. IIVI 


II. BACKGROUND 


For simplicity, we assume that only the highest occu¬ 
pied molecular orbital (HOMO) contributes to the dy¬ 
namics. We also neglect the motion of the nuclei and 
represent the HOMO by a linear combination of atomic 
orbitals (LCAO). In this case, the HOMO wavefunction 
'Fo(r) = (r | H/o) is written as 


^o(r) =^< 


Ipa 







(3) 

where ip a ( r) are the atomic orbitals, R is the internuclear 
distance, £ a is the orbital quantum number and m a is the 
magnetic quantum number. The indices A a = m a relate 
to gerade (g) and A a = m a + 1 to ungerade (u) orbital 
symmetry, respectively. We incorporate the molecular 
structure in the prefactors and keep the action unaltered. 
This is a reasonable assumption if the internuclear dis¬ 
tance is not very large (see, e.g., (' [23 - |26| ) for details). 

The transition amplitude © is computed using the 
steepest descent method, i.e., we seek t', t and p so 
that the action ED is stationary. These stationary so¬ 
lutions can be associated with the tunneling time t' and 
the rescattering time t of an electron propagating in the 
continuum, and it constrains its intermediate momentum 
p so that it returns to its parent ion [ 7 ]. We employ either 
the uniform approximation or the saddle-point approxi¬ 
mation discussed in Ref. [27 ]. 

We will now assume that the field is made up of two 
linearly polarized orthogonal laser fields. Explicitly, the 
time dependent electric field E(t) = —dA(t)/dt and the 
vector potential A (t) may be written as 


A. Single-molecule response and orbit-dependent 
shift 


The SFA transition amplitude for HHG [t| reads as 

/ OO pt /* 

dt dt' / d 3 pd* ec {p + A(t))d io n(p + A(f')) 
- 00 J — OO J 

xe iS(t,t',n, P ) + C-Cj (!) 


where d ion (p) = (p|R/(i')l^o> and d rec (p) = (pM’I'o) 
are the ionization and recombination dipole matrix ele¬ 
ments, respectively, and 

1 /■* 

S{t,t',fl,p) = —- J' [p + A(T)] 2 dT-I p (t-t') + £U (2) 


denotes the semiclassical action. Eq. ED describes a pro¬ 
cess in which an electron propagates in the continuum 
from the ionization time t' to the recombination time t 
with field-dressed momentum p + A (r) , i' < r < t. 
The ionization potential, the vector potential and the 
harmonic frequency are denoted by I p , A and H, respec¬ 
tively. Our computations have been performed in the 
length gauge, so that the interaction Hamiltonian is given 
byF;(0 = r-E(t'). 


E(f) — E\\ (f)ey + Ej_(t)e± 

(4) 

A (t) = A||(t)e|| + Aj_(f)e_L, 

(5) 


where the field components and unit vectors along the 
major and the minor polarization axis are indicated by 
the subscripts (||) and (_L), respectively. 

The structural interference patterns are mainly given 
by the recombination prefactor. Here, we are interested 
in d* ec (p(t) ■ E(t)), which is the complex conjugate of its 
component along the field-polarization direction. Explic¬ 
itly, 

d* rec (p(t) • E(t)) = ^c a [e- ip W'^ + (-i)G— «+A 0e ip(t)-f 

a 

x(-f)^I a pnbV’a( p{t)) E b(t), (6) 

b 

where p (t) = p + A(£), 

= ( 2?r ) 3/2 J d 3 rip a (r)exp[-ir-p(t)}, (7) 
and b= 11, J_ 
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After lengthy, but straightforward computations, one 
may show that interference minima will occur for har¬ 
monic order 


n =_ 2[nn ~ a]2 _+ / 

R? CO s 2 (0 L - C(M')) P 


( 8 ) 


where 6l is the alignment angle of the molecule with 
regard to the main polarization axis, and 


((t, t') = arctan 


P. l + A±(t) ' 
.P\\ + A \\(.*) \ ’ 


(9) 


is an effective shift caused by the non-vanishing field el- 
lipticity. In Eq. ©, 


Pb 


-1 

t-t' 



A b (r)d,T , 


( 10 ) 


with b =|| and b =_L, denote the intermediate momentum 
of the electron along the major and minor polarization 
axis, respectively. 

In other words, the angle of return has effectively 
been incorporated in the two-center interference condi¬ 
tion. The shift a = arctan(— iA+/A-) is structural, and 
related to s — p mixing, with 

-4±=E c “ ± 1] b(P + A(i), t) (11) 

a 

and 

V(P,t) = ~i [d Pl] ip* a (p)E n (t) + d p± il>Z{p)E ± (t)] . (12) 

For details on this derivation we refer to our previous 
publication Q For linearly polarized fields, £(£, t') = 0 
and Eq. © reduces to the purely structural condition 
derived in 1281 . 


B. Propagation model 


To model experimental data, we must consider the 
macroscopic harmonic response generated by an intense 
laser pulse focused into a gaseous medium, typically sup¬ 
plied by a gas jet directed perpendicularly across the 
beam. We must then numerically integrate Maxwell’s 
wave equations with source terms distributed across this 
extended medium. After applying the slowly evolving 
wave approximation (SEWA) d 2 E/dz 2 = 0, Maxwell’s 
equations can be expressed in frequency space as [29] 


d . ic _ 9 . 2niuj -. . 

— £l(r, z, w) + —V 2 ± E(r, z, w) =- P(r, z, w), 

oz 2cj c 

(13) 

where E(r,z,w) is the electric field comprised of both the 


driving laser field and the generated harmonic field, and 
P(r, z , w) is the polarisation response of the medium. We 


assume that the radiation generated by nonlinear inter¬ 


actions does not influence the strong driving field, which 


is a good approximation for the parameter range of inter¬ 
est. This means that we can separate the total electric 
field into two distinct components: the driving IR laser 
field Ei(r,z,oj) and the generated XUV field Eh{r,z,uj), 
allowing us to solve Eq. m separately for each compo¬ 
nent using the appropriate approximations for each case. 
We consider atomic gas densities between 10 16 and 10 17 
atoms/cm 3 and interaction lengths of approximately a 
few mm. Under these conditions we can ignore all linear 
dispersion and absorption effects for the driving field and 
consider only the polarization response Pi on (r, z,w) due 
to the oscillation of the free electrons created through 
ionization of the gas medium. 

For the XUV component of the propagating field 
we can ignore free electron effects as these frequencies 
quickly exceed the plasma frequency. However, we must 
include the nonlinear dipole response of the molecular 
gas P(r,z,uj) and XUV absorption by the neutral atoms 
via the frequency dependent absorption coefficients a a bs 
of Ref. [30]. In the present calculations linear dispersion 
due to neutral atoms has not been included. This sim¬ 
plification is justified for H 2 , as this molecule exhibits no 
resonances in the continuum for the frequency range of 
interest. 

The contribution from the nonlinear dipole response is 
advanced according to 

■^~E h (r,z,uj) = -^^-P(r,z,uj), (14) 

oz c 

where 


P(r,z,u)) = n a {z)Mz(r,z,u>). (15) 

In Eq. (ED, Mi(r,z,bo) is the frequency spectrum of the 
dipole acceleration along e, calculated using the single 
atom response model of Sec. Ill Al and n a (z) is the gas 
density, which we have defined using the Gaussian profile 

n a (z) = N a e~ 4ln2( '^ L)2 . (16) 

In the above-stated equation, N a is the peak gas density, 
Zj marks the z position of the centre of the gas jet, and 
z u is the full width at half maximum (FWHM) width of 
the jet. 

The far-field response, which is actually observed in ex¬ 
periments, may be obtained analytically in the frequency 
domain via a Huygen’s integral. This gives 

Eh(r, z + lz,u>) =^- [ r , E h (r,z,u})e~^^ 2+r ' 2 ' ) 

« J 00 

VV f UJ 

J ° ( ^f )dr '’ (17) 

where l is the distance of propagation to the far field 
and Jo is there zeroth-order Bessel function. For more 
details on the propagation model we refer to our previous 
publication @]. 

The measured XUV field is the coherent sum of all the 
harmonics generated from all the atoms in the nonlinear 
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medium. The far-held spectrum is composed of those 
components of the single atom dipole emission, from spe¬ 
cific regions of the focal volume, for which there exists 
a wave vector along which the components add coher¬ 
ently [ 20 !, Hi . The degree of phase matching is usually 
dominated by the spatial variation of the intensity and 
carrier-envelope phase (CEP) of the laser held through 
the interaction region. Therefore, changing the gas jet 
position relative to the laser focus or changing the gas jet 
density prohle can drastically change the phase-matching 
conditions. Most notably, it can be used to select either 
short or long plateau contributions [HI HI or the half¬ 
cycle cut-off emissions [23|. 


III. RESULTS 



In the results that follow, we consider H 2 in a held 
composed of the orthogonally polarized waves 

E\\ = E u f{t) cos(wi) (18) 


and 


E _l = E 2 uif(t) cos(2wf + (/)), (19) 

where the pulse shape f(t) is a Gaussian function. Unless 
otherwise stated, the fundamental pulse ( 1181 ) has wave¬ 
length A = 800nm, peak intensity I u = 2.5 x 10 14 W/cm 2 , 
and a full width at half maximum of 30 fs (approximately 
10 cycles), which are within the experimentally relevant 
parameter range. 

For H 2) the HOMO is a \<r g orbital composed of s 
orbitals only. Hence, the interference condition reduces 
to 


n = i p 


2n 2 TT 2 
R 2 cos 2 0l 


( 20 ) 


for linearly polarized helds. For nonvanishing ellipticity, 
6l must be replaced by 9l — in the above-stated 

equation. The orbitals employed in this work have been 
computed using GAMESS-UK jH}. 


A. Intensity dependence 

We will hrst get an insight on how the phase shift 
£(t, t') depends on the driving-field intensity. Since this 
parameter varies strongly across the laser beam profile, a 
simplified preliminary investigation is useful in order to 
understand its overall behavior. To facilitate the inter¬ 
pretation, we will approximate the pulses ( 1181 ) and m 
by monochromatic waves. 

The shifts depend strongly on the ionization and re¬ 
combination times, which vary across the spectra. Specif¬ 
ically, there is a dominant pair of orbits, the well-known 
“long orbit” and “short orbit”. For the long orbit, the 
electron leaves slightly after the field peak and returns 


FIG. 1: (Color online) Real parts of the effective shifts ( ( t , t') 
as functions of the driving-field intensity of the fundamen¬ 
tal for a harmonic of frequency Q. = 31a;, using H 2 (J p = 0.59 
a.u.) in two-color laser fields of increasing ellipticity. The fun¬ 
damental and the second harmonic have been approximated 
by monochromatic waves where bu = The ellipticities 
have been increased from £ = 0 to £ = 0.3 in increments of 
= 0.05. A lighter color indicates a higher ellipticity. The 
dashed and solid lines refer to the short and long orbits, re¬ 
spectively. Panel (a) and (b) refer to relative phases of <f> = 0 
and 4> = 0.4-7T, respectively. 


after a field crossing, while for the short orbit it ionizes 
at a later time and returns before the crossing. The lower 
the harmonic is in the plateau, the further apart these 
times are. In the cutoff region, these orbits coalesce. 

In Fig. [L] the real parts of are plotted for a 

specific harmonic, as functions of the driving-field inten¬ 
sity. The figure shows that the shifts depend strongly 
on the held ellipticity, on the phase difference cf> between 
the two waves and on the driving-held intensity. This 
is expected as variations in the intensity will change the 
energy position occupied by a specihc harmonic. This is 
related to the fact, that, for low enough intensities, these 
shifts nearly coincide for the long and short orbits. In 
this case, the harmonic employed is located either at or 
beyond the cutoff frequency. As the intensity increases, 
there is a splitting in the shift, indicating that the plateau 
region has been reached. 

The behavior with regard to the field parameters de¬ 
pends very strongly on the orbit. If the two waves are in 
phase, for instance [Fig. Oja)], the distinguishing features 
are a large residual shift after the cutoff and a very pro¬ 
nounced negative shift for the short orbit in the plateau 
region. Physically, this means that if the electron is re¬ 
turning along the short orbit, for <fi = 0 the angle of 
return is much larger. Hence, it may be easier to single 
out in a realistic scenario, especially at high intensities. 

This behavior, however, changes with the relative 
phase. An example is provided in Fig. [jjb), which ex- 
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hibits a smaller residual shift in the cutoff region and 
comparable shifts of opposite signal if the chosen har¬ 
monic is in the high-plateau energy region. Increasing 
the intensity will lead to a non-trivial behavior for the 
shift associated with the short orbit, while the shift re¬ 
lated to the long orbit will exhibit a monotonic behavior. 
Physically, this implies that the electron’s angle of return 
will vary much more dramatically across the beam for the 
short orbit. 


B. Spatial effects across the beam profile 


We will first establish how the structural minimum 
manifests itself in the macroscopic case, and then dis¬ 
cuss the situation with fields of non-vanishing ellipticity. 
Throughout, we have chosen the propagation conditions 
so that the short orbit is favored. In general, this is 
achieved by placing the center of the gas jet after the fo¬ 
cus. This type of configuration leads to an enhancement 
of on-axis harmonic emission, with low divergence. Fur¬ 
thermore, we have chosen the ionization times to start 
shortly after the central maximum of the pulse, so that 
the corresponding return times will occur near the sub¬ 
sequent crossing. 

In Fig. [2j we examine the structural minimum across 
the beam for a linearly polarized field, both near the in¬ 
teraction region and in the far field. This minimum is 
indicated by the straight vertical lines, which correspond 
to the interference condition Eq. (l20l) and agree well with 
the suppression in the yield. The energy position re¬ 
lated to this suppression remains the same regardless of 
whether contributions of individual orbits are taken, as 
shown in panels (a) and (b), or if they are combined, as 
shown in the remaining panels. Furthermore, it does not 
change with the driving-field intensity [panel (c)] or in 
the far-field regime [panels(e) and (f)]. Only when the 
alignment angle is varied, namely for panels (d) and (f) 
does this minimum change. All this is consistent with the 
fact that the two-center minimum described by Eq. (l20l) 
is purely structural. The fringes in panels (c) to (f) are 
caused by the interference of the long and short orbits. 

In the upper panels of Fig. [3] we present the macro¬ 
scopic HHG spectrum computed for orthogonally polar¬ 
ized fields considering a coherent superposition of the 
long and short orbits, and a phase difference of </> = — 0.27T 
between the to and 2w waves. This figure illustrates a par¬ 
ticular case for which the shift associated with the short 
orbit, given by the red lines, is visible in a realistic sce¬ 
nario, and survives in the far field. The close-up near 
the interference minimum [Fig. (3][b)] shows that it fol¬ 
lows the generalized interference condition related to the 
short orbit very closely. Only at the beam edges there 
are small discrepancies from this condition, and the sup¬ 
pression approaches the shifted minimum associated to 
the long orbit, which is given by the green line. 

The reason behind this clear picture can be seen in 
panels (c) and (d) of Fig. [3l in which we present the 


FIG. 2: HHG macroscopic response of H 2 in a linearly polar¬ 
ized field (left and right panels, respectively) of wavelength 
A = 800nm, plotted in a logarithmic scale. The beam waist 
is w = 30/rm, and the gas jet is placed at z g = 4mm after the 
focus. The FWHM of the gas jet is 0.5mm and the FWHM of 
the intensity envelope is 30fs. Panels (a) and (b): Individual 
contributions of the long and short orbits, for a driving-field 
intensity /u, = 2 x 10 14 W/cm 2 , and alignment angle 9 = 0 ; 
note that the oscillations in Panel (b) are caused by an artifact 
due the inaccuracy of the standard saddle-point approxima¬ 
tion in the cutoff region. Panels (c) and (d): spectra from 
the coherent superpositions of the long and short orbits, for a 
driving-field intensity I^ = 2.5 x 10 14 W/cm 2 , and alignment 
angles 9 = 0 [panel (c)] and 9 = 7 t/6 [panel (d)[; panels (e) 
and (f): far-field harmonic spectra, for the same intensity and 
alignment angles as in panels (c) and (d). The dashed line 
indicates the position of the structural two-center minimum. 
The labels s, l and l + s indicate contributions from the short 
orbit, long orbit, or from a coherent superposition of both, 
respectively. 


individual contributions from the long and short orbits. 
They show that the contributions from the long orbit are 
strongly suppressed in the plateau region, and that those 
from the short orbit are much more significant. This 
stems from the fact that the relative phase </> = — 0.27T 
selects the short orbit at the single-molecule level via po¬ 
larization gating Q • This selection is then reinforced by 
the appropriate propagation conditions. Also for that 
reason there are no visible fringes when the coherent su¬ 
perposition of both orbits is considered (Fig. [3]). 

Next, we will analyze the generalized interference con¬ 
dition for individual orbits in more detail. With that 
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FIG. 3: Propagated HHG spectra for H 2 in a Gaussian pulse, 
where the peak intensities of the uj and 2<u waves are I u = 
2.5 x 10 14 W/cm 2 and l 2 U = 7.5 x 10 13 W/cm 2 defined at 
the gas jet (their intensity ratio is around 0.3) and relative 
phase <j> = — 0.2-7T. The wavelength of the fundamental is A = 
800 nm. The beam waist is w = 30 fim and the center of the 
gas jet is located at z g = 4 mm after the focus. The FWHM of 
the intensity envelope is 30fs. Panels (a) and (b) display the 
far-held result for a coherent superpostion of orbits denoted 
by the labels l + s, where (b) is a close-up of the region around 
the shifted minimum seen in (a). Panels (c) and (d) show the 
individual contributions to far held spectra for the long and 
short orbit respectively. The dashed lines and the solid lines 
indicate the position of the unshifted and shifted two-center 
minima, respectively. The red and the green solid lines give 
the positions of the shifted minima for the short and long 
orbits, respectively. The curves have been computed spread 
equally across the interaction region such that z g — Az < z < 
z g + Az, with Az = 0.5 mm. All spectra have been plotted in 
arbitrary units and in a logarithmic scale. 


aim in mind, we (i) choose a phase difference 0 = 0, for 
which these shifts are expected to be large EH ; (ii) con¬ 
sider the near-field regime so that diffraction effects are 
ruled out. These results are plotted in Fig. |1J for the 
short and long orbits. The figure shows that these shifts 
depend very strongly on the orbit. While, for the long 
orbit, the shift moves the minimum away from the cut¬ 
off and towards lower frequencies, for the short orbit it 
causes it to move towards the cutoff [see panels (a) and 
(b)]. Furthermore, a zoom-in of the shifts, presented in 
panels (c) and (d), shows that they approach each other 
around the cutoff and follow the modified interference 
condition f5]) very closely. In particular, the shift is much 
larger for the short orbit and, as long as the cutoff has 
not been reached, it varies very little within the interac¬ 
tion region. For the long orbit, the shift varies slightly 
more with regard to z. Beyond the cutoff, the variation 
is more extreme. This region, however, is not of interest 
to the present problem. 

Unfortunately, however, these shifts cannot be seen if 


FIG. 4: Individual contributions of the long and short orbits 
(left and right panels, respectively) for the HHG macroscopic 
response of H 2 in a Gaussian pulse composed of parallel and 
perpendicular waves given by Eq. m and m, respectively. 
The peak intensities of the uj and 2ui waves are I u = 2.5 x 
10 14 W/cm 2 and hu, = 7.5 x 10 13 W/cm 2 defined at the gas jet 
(their intensity ratio is around 0.3) and relative phase <j> = 0. 
The center of the gas jet is located at z g = 2 mm after the 
focus. This phase has been found to give a very large residual 
shift ((t, t') at the cutoff. The wavelength of the fundamental 
is A = 800nm. The beam waist is w = 30[im, and the FWHM 
of the intensity envelope is 30fs. Panels (a) and (b): HHG 
yield from the individual orbits; panels (c) and (d): zoom in 
of the upper panels close to the interference minimum. The 
dashed lines and the solid lines in the upper and middle panels 
indicate the position of the unshifted and shifted two-center 
minima, respectively. The red and the green solid lines give 
the positions of the shifted minima for the short and long 
orbits, respectively. The curves have been computed spread 
equally across the interaction region such that z g — Az < z < 
z g + Az, with Az = 1 mm. All panels have been plotted in 
arbitrary units and in a logarithmic scale. The labels l and s 
are associated to the long and the short orbit, respectively. 


a coherent superposition of orbits is considered. We have 
plotted this superposition in Fig. O for the near and far 
field regimes [upper and lower panels in the figure, re¬ 
spectively]. The overall behavior shows many interfer¬ 
ence fringes and no clear shifted minimum. This happens 
because, at the single-molecule response level, a phase 
0 = 0 will enhance the long orbit and suppress the short 
orbit Q. Hence, there are conflicting conditions from 
the the single-molecule response and propagation, which 
is not ideal. This is clearly seen in panels (a) and (c) 
of the figure, which exhibit fringes in the whole plateau 
region due to the interference between the short and long 
orbits. Close to the two-center minimum the results are 
inconclusive as (i) the suppression near the green curve 
seems much more related to the interference between the 
long and short orbit than to the electron’s angle of return 
being incorporated; (ii) the shift associated to the short 
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FIG. 5: Propagated spectra considering a coherent superposi¬ 
tion of the long and short orbits, for the same field parameters 
as in Fig. [4] Most propagation conditions have also been kept 
as in Fig. [I] except the center of the gas jet, which has been 
chosen to be at z g = 2 mm and z g = 4 mm after the focus 
(left and right panels, respectively). Panels (a) and (b) dis¬ 
play the spectra in the interaction region, while panels (c) 
and (d) show the far-held results. The dashed lines and the 
solid lines indicate the position of the unshifted and shifted 
two-center minima, respectively. The red and the green solid 
lines give the positions of the shifted minima for the short 
and long orbits, respectively. The curves have been com¬ 
puted spread equally across the interaction region such that 
z g — Az < z < z g + Az, with Az = 1 mm. All spectra have 
been plotted in arbitrary units, and in a logarithmic scale. 

orbit does not manifest itself as a clear suppression. Both 
issues (i) and (ii) become slightly better in the far-held, 
as the interference between the short and long orbit is 
partially washed out. Nonetheless, the result is not as 
clear as in Fig. [3j for which the short orbit was favored 
at the microscopic and macroscopic level. 

This problem may be attenuated by moving the gas jet 
further away from the focus, in order to favor the short 
orbit [see Figs. (5](b) and (d)]. In this case, the interference 
fringes in the plateau become more blurred, so that the 
shift related to the short orbit can be identified. 

One should bear in mind, however, that we have re¬ 
stricted the contributing orbits to a single half cycle of the 
driving field. This does not correspond to a realistic situ¬ 
ation, as a real pulse would have many cycles. A specific 
problem is that, in a bichromatic u>—2u> field with orthog¬ 
onally polarized waves, the long and short orbits starting 
at the subsequent half cycle gives rise to shifts of 

opposite signs 0 - In practice this causes a blurring in 
the generalized interference condition. This problem may 
be overcome if one employs a few-cycle pulse. 

This is shown in Fig.[6l in which we present the far-held 
spectra calculated using the two dominant half cycles of 
the pulse in Fig. [3l and of a few-cycle pulse with the 
same parameters, except the full width at half maximum 


FIG. 6: Propagated HHG spectra considering two half cycles, 
for the same driving-field parameters in Fig. [3] but different 
FWHM of the intensity envelope. In panels (a) and (c), this 
width is 30fs, while in panels (b) and (d) it is 5.5fs. In the 
lower panels we show only the shift related to the short orbit, 
which dominates throughout. These shifts are displayed as 
the thick lines in the figure. The shifts associated with the 
first and second half cycle are shown as the red curves in pan¬ 
els (c) and (d), and the pink curve in panel (c), respectively. 

of the intensity envelope (left and right panels, respec¬ 
tively). For the long pulse, the contributions from both 
half cycles blur the minimum associated with the short 
orbit [panel (c)], while for the few-cycle pulse it remains 
very clear. It is also well described by the interference 
condition © [panel (d)]. This happens because the two- 
center minimum is in the plateau for the first half cycle, 
while it is beyond the cutoff for the second half cycle 
of the short pulse. The different cutoff energies become 
clear in panel (b), which shows well-defined harmonics 
only for energies lower than the structural minimum. In 
contrast, the long pulse leads to well-defined harmonics 
throughout, as seen in panel (a). 


IV. CONCLUSIONS 

In this work, we have performed theoretical studies 
of the macroscopic response of high-harmonic radiation 
from diatomic molecules in bichromatic fields composed 
of two orthogonally polarized driving waves. Our com¬ 
putations show that, for an appropriate choice of driving 
field and propagation conditions, the angle of return of 
an electron to its parent molecule manifests itself as a dy¬ 
namic shift in the two-center structural minimum. This 
shift depends on the driving-field intensity, on the har¬ 
monic order and on the orbit along which the electron 
returns. For the dominant orbits, these shifts are visible 
both in the near- and far-field regimes, indicating that 
they can, in principle, be measured in experiments. 






























In order to see these shifts clearly, one must select an 
individual orbit by employing gating mechanisms at the 
single molecule level and by using appropriate propaga¬ 
tion conditions. It is of particular interest to enhance 
the short orbit and suppress the long orbit, as it phase 
matches on axis and leads to more dramatic variations 
in the angle of return with the laser-field parameters. If 
conflicting microscopic and macroscopic conditions are 
provided, the interference between long and short orbits 
may lead to inconclusive results. With that aim in mind, 
we have selected the phase difference between both ui and 
2w driving waves as <f> = — 0.27T and have placed the focus 
before the gas jet. These choices are known to enhance 
the short orbit at the single-atom [§J and macroscopic 
mill levels, respectively. 

Nonetheless, there exist two obstacles towards seeing 
this shift, which, however, are well under control. First, 
for bichromatic fields, the shifts coming from the other 
half cycle will have opposite sign. This flip may be 
avoided for few-cycle pulses, by creating a spectral and 
temporal gate in order to select electron orbits from a 
single half cycle [23|, [34j ]. This can be performed if for 
one of the dominant half cycles, the shifted minimum 
lies beyond the cutoff energy, and within the plateau for 
the other half cycle, as shown in Fig. [6] 

Second, for molecules with a more complex electronic 
structure, the core dynamics may also play an important 
role, such that the single-active orbital, single-active elec¬ 
tron approximation is not applicable [161 1 • The regime for 
which this happens, however, can be avoided. Indeed, 
recent experiments in CO 2 have shown that interference 
effects stemming from the core dynamics are only rele¬ 
vant if the structural minimum provides a window for 
them to be observed. For high enough intensities, the 
energy position of the dynamical minimum lies outside 
this window and can no longer be observed [3f|. This 


is in agreement with the experimental findings in [36j |. 
which support the structural instead of the dynamical 
interference picture, and with the computations in [37} . 
which show that multielectron effects are important in 
the cutoff, but not in the mid-plateau energy region. 

Finally, our results suggest that, in principle, the an¬ 
gle with which the electron returns to the core can be 
controlled, either by changing the relative phase between 
both waves, or the driving-field intensity. As the inten¬ 
sity increases, a specific harmonic moves from the cutoff 
region across the plateau towards the ionization thresh¬ 
old, and this will alter the angle of return. This behav¬ 
ior is particularly critical for the short orbit, which, as 
discussed above, is the most favorable for macroscopic 
observations of this shift. Furthermore, this angle can be 
mapped into a shift near the structural minimum, which 
can be modified by an appropriate parameter choice. 
This may be desirable in future experiments in order to 
investigate dynamic effects which would be obfuscated 
otherwise. 
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